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ABSTRACT 

We have developed a perturbative approach to microlensing due to an 
extrasolar planetary lens. In particular, we have found analytic formulae 
for triple images. We have used the formulae to investigate the astrometric 
microlensing due to the extrasolar planetary lens, in expectation of dramatic 
improvements in the precision of the future astrometric measurements. For 
a weak lensing case, we have shown how the maximum angular size and the 
typical time scale of the anomalous shift of the light centroid are dependent 
on the mass ratio and angular separation between the star and the planet. 

Subject headings: gravitational lensing — astrometry — planets 
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1. Introduction 

Extrasolar planets searches are successfully going on (Marcy & Butler 1998). The 
number of candidates for the planets has reached about 70 thanks to the Doppler 
method. In spite of the success, this technique has some limitations: The inferred mass 
is the lower bound, since the inclination angle of the orbital plane is not determined 
except for the echpse case. In addition, the radial velocity of a star becomes too small to 
detect its Doppler effect when separation between the star and the planet is of the order 
of an AU. Hence, it seems quite difficult to discover earth-type planets by this method. 
Astrometry is considered as a supplementary method to find out planets by measuring 
the transverse motion of a star. Indeed, we can expect significant improvements in the 
precision of the future astrometric measurements by DIVA, FAME, GAIA and SIM. 
Their precision will achieve a few micro arcseconds. 

These future astrometry missions make the astrometric microlensing the third 
method for searching extrasolar planets: The typical scale in the gravitational lensing 
is the Einstein ring radius, which is of the order of an AU for a lensing star within our 
galaxy. Hence, the microlensing is quite effective even for earth-type planets (Mao & 
Paczynski 1991, Gould & Loeb 1992). In the planetary lens case, we can determine its 
true mass through observation of spikes in the light curve, which are produced by the 
magnification effect of the gravitational lensing. On the other hand, the astrometric 
microlensing is a consequence of a combination of the magnification and the position 
shift of the image (Miyamoto and Yoshii 1995). The numerical results (Safizadeh et al. 
1999, Han and Lee 2002) have shown that the photo centroid shifts provide us a clue 
for extrasolar planets. However, the parameter dependence is not clear in the numerical 
approach. The purpose of this paper is to derive the analytic formulae and clarify the 
parameter dependence. Analytic approaches to a binary lens have not been developed 
until now. After presenting a brief summary of the astrometric microlensing due to a 
single lens, we treat a planetary lensing case as a perturbation around the stellar lens. 
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2. Microlensing due to a planeteiry system 

2.1. A single lens 

A single lens is located at the distance Dl from the observer, and a source at Ds- 
The distance between the lens and the source is denoted by Dls- Under the thin lens 
approximation, the lens equation for the single lens is written as 

/3 = 9-^a, (1) 

where /3 and are the angular positions of the source and the image, respectively, and oc 
is the deflection angle. For a spherical lens with mass M, or at least axisymmetric along 
the hne of sight, the deflection angle is 

Ia| = (2) 

where G is the gravitational constant and b is the impact parameter of the light ray. It 
is convenient to normalize our equations by the angular radius of the Einstein ring 



_ 4GM^ 



which is of the order of a milli arcsecond (mas) for a solar mass lens within our galaxy. 
In units of the Einstein ring radius, the lens equation is rewritten as 

/3 = ^-^, (4) 
where 6 denotes \6\ and similar notations are used below. We find out the two solutions 

= A;(±)/3, (5) 

where we defined 

The amplification for each image is respectively 
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Following the definition of the center of mass, we define the photo center as (Walker 
1995) 

^'^ AM + M-) • 

With respect to the unlensed position, the location of the photo center is 

Mc = 0c~f3, (9) 

which expresses the deviation due to the lensing. For the single lens, it becomes simply 

A«c = ^. (10) 

where we used Eq. (0). We can assume that the stellar motion in our galaxy is 
approximated over several decades by a straight line, since the curvature of the orbit in 
our galaxy is negligible. For the later convenience, we decompose f3 into the orthogonal 
vectors as 

(3 = b + V, (11) 
where b is the vector for the impact parameter from the lens to the stellar orbit, and 

V = tv^, (12) 

for Vj_, the transverse angular velocity of the source to the lens. Here, we have chosen 
t = as the time when the star is closest to the lens. We introduce C, as 



V = ^JP^ + 2tan^, (13) 

so that the photo center can be rewritten as 

b bcos2e Vsin2e 

^ 2(62 + 2) 2{b^ + 2)^ 2VV^T2' ^ ^ 

which is an ellipse (Walker 1995, Jeong et al. 1999). Its size is of the order of a mas 
for stellar cases in our galaxy. The center of the ellipse is at 2{b^+2) ' length of 

the semimajor and semiminor axes are 2^b'^+2 2(b^+2) ' respectively. The photo center 
passes the two points on the semiminor axes at t = (^ = 0) and t = ±oo (^ = ±f ). 
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2.2. A planetary lens 

Now, we are in the position to consider the astrometric microlensing due to the 
planetary system, where the mass of the star and the planet are Mi and M2 respectively. 
Let us consider the separation vector from the star to the planet. Its projection onto 
the lens plane is denoted by s. We adopt the frame of center of mass. Then, in the unit 
of the Einstein ring radius angle due to the total mass Mi + M2, the lens equation is 
written as 



where we defined 



e - (18) 

In the planetary case, M2 is much smaller than Mi. For instance, the Jupiter mass 
is about 10~^ of the solar mass. Hence, we introduce an expansion parameter as u — 1^2 
in our perturbation approach. Since we wish to consider the microlensing as a method 
supplementary to the Doppler technique, we concentrate on a large separation case 
e > 1, which is beyond the reach of the Doppler method. It is straightforward to extend 
our investigation to the case of e < 1. In addition, we consider a case of a large impact 
parameter /5 3> 1, which is most probable because of the large cross section. In total, 
we consider the case of /3 ^ e > 1. Let us look for the solutions of the lens equation by 
taking a form of 

6^60 + 66, (19) 
where Oq and SO are the solutions at the zeroth and first order, respectively. 
We know the zeroth-order solution 

0f ^ = /c(±)/3, (20) 
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which is expanded as 

(21) 

Ot^ ~ « e. (22) 

That is, the "+" image appears outside the planetary system, while the "— " image 
inside it. Throughout our calculations, we should keep in mind that the large impact 
parameter does not mean 9 ^ e. 

Next, we find the solutions at the first order as 

<Se<+' = ^ (-2e(/3 . e)0' - fie'ffi + 4/3(/3 • e)^ + 0{^l3^)) , (23) 
= , (_e+^ + 0(-L)), (24, 

where the first and second solutions are perturbations around the zeroth-order solutions, 
and the third solution does not appear until at the linear order of u. Actually, for z/ = 0, 
the third image is located at the direction to the planet. As for the number of the 
images, see also Witt (1993). 

Hence, we obtain the Jacobian of the mapping between the source and the lens 
planes. 



,1^6 



4 



= l + 0(-^,^^), (26) 



' a/3 ' 

1^1 = -— + 0(.^), (27) 
I^HO(.^). (28) 

The contribution of the third solution to the photo center is negligible for the large 
impact parameter. 
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2.3. Distortion of an ellipse for light centroid shifts 

By the use of the preceding subsection, we find a correction to the photo centroid as 

5dc = ;^(-26(/3 ■ e)(3' - f3e'(3' + 4(3i(3 ■ ef + 0{(3')) + 0{v^). (29) 

This is due to a primary effect caused by a position shift given by Eq. (p3D, since Eq. 
(p6|) shows that the correction to amplification appears only at higher orders. 

Let us investigate in detail the photo center shift. We consider the expansion of the 
following superposition of two ellipses 

13 + TJe (3-e 

(1 - , ..19 , o + ^ 



|/3 + //e|2 + 2 |/3-e|2 + 2 
/3 1/ 



+ -5^(-2e(/3 ■ e)(3^ - f3e^(3^ + 4/3(/3 ■ e)^ + 0{e'(3^)) + 0{u^), (30) 



/?2 + 2 /36 

where the first term of the left hand side is considered as the primary ellipse and the 
second as the perturbation around the primary. The position of the first ellipse is shifted 
by ve in comparison with Eq. (|lO|) , and its size changes by factor 1 — v. Next, let us 
take a closer look at the secondary ellipse. The orthogonal decomposition by Eq. ([TID is 
modified as 

/3 - e = (6 - ex) + (t - tc)^±, (31) 

where we defined 

^„ - (32) 

e_L = e — e|| , (33) 

tc = ^. (34) 

For stellar cases in our galaxy, the maximum angular size of the distortion is 
estimated as 

^(i)' ^ micro arcsec.(^)(i)(^)(^)^ (35) 
where we assumed that (3 is comparable to e. The maximal distortion occurs at 

,_10M^)(^). (36) 
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about a few months before/after, depending on a location of the planet, the source passes 
the point closest to the lensing star. A planet whose orbital separation is between 0.1 
AU and a dozen of AU will be detectable by future missions which observe at intervals 
from a few days to several years. 



2.4. Discussion 

The position shift falls off as b~^, while the amphfication as Hence, the cross 
section for the astrometric microlensing is much larger than that for the photometric 
microlensing, when the sufficient accuracy of astrometric measurements will be achieved. 
This is why, up to this point, we have concentrated ourselves within a case of a large 
impact parameter. In a small impact parameter case, the caustic crossing produces 
sufficiently large magnification of the images, so that we might be able to detect it 
much more easily by the astrometric microlensing as well as photometric microlensing 
(Safizadeh et al. 1999). Numerical implementations are needed to study such a nonlinear 
behavior. In order to evaluate feasibility in the future mission, we must take account 
of the brightness of the lensing star. It might be important to pay attention also to 
fluctuation of extragalactic reference frame due to gravitational lensing of black holes 
(Schutz 1982) and MACHOs (Hosokawa et al 1997) in our galaxy. Finally, the Keplerian 
motion might cause an appreciable effect on the distortion of the photo center ellipse, 
since the orbital period of the planet ranges from months to more than years, presumably 
comparable to to- Therefore, it would be an important subject to study these effects in 
detail. 



3. Conclusion 

In expectation of dramatic improvements in the precision of the future astrometric 
measurements, we have developed a perturbative approach to microlensing due to an 
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extrasolar planetary lens. Formulae for triple images due to the planetary lens are given 
perturbatively by Eqs. (p3D — (|25|). In particular, we have shown by Eqs. (|35|) and (|36D, 
how the light centroid shifts are dependent on the mass ratio and separation between the 
star and the planet. The typical time scale is of the order of months, depending strongly 
on ey, a projection of the separation vector onto the source trajectory. 
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